Draft version December 18, 2012 

Preprint typeset using I^TjrX style emulateapj v. 5/2/11 



GAS AND DUST ABSORPTION IN THE DOAR 24E SYSTEM 
Andrew J. Kruger 

Department of Physical Science, Wilbur Wright College, 4300 N. Narragansett Ave, Chicago, IL 60634 



Matthew J. Richter 

Department of Physics, University of California at Davis, One Shields Ave, Davis, CA 95616, USA 

Andreas Seifahrt 

Department of Astronomy and Astrophysics, University of Chicago, 5640 S. Ellis Ave, Chicago, IL 60637, USA 

John S. Carr 

Remote Sensing Division, Naval Research Laboratory, Code 7210, Washington, DC 20375, USA 

Joan R. Najita 

National Optical Astronomy Observatory, Tucson, AZ 85719, USA 

Margaret M. Moerchen 

European Southern Observatory, Alonso de Cordova 3107, Santiago, Chile and 
Leiden Observatory, PO Box 9513, 2300 RA Leiden, The Netherlands 

AND 

Greg W. Doppmann 

W.M. Keck Observatory, 65-1120 Mamalahoa Hwy, Kamuela HI 96743 
Draft version December 18, 2012 

Abstract 

We present findings for DoAr 24E, a binary system that includes a classical infrared companion. We 
observed the DoAr 24E system with the Spitzer Infrared Spectrograph (IRS), with high-resolution, 
near-infrared spectroscopy of CO vibrational transitions, and with mid-infrared imaging. The source 
of high extinction toward infrared companions has been an item of continuing interest. Here we 
investigate the disk structure of DoAr 24E using the column densities, temperature, and velocity 
profiles of two CO absorption features seen toward DoAr 24Eb. We model the SEDs found using 
T-ReCS imaging, and investigate the likely sources of extinction toward DoAr 24Eb. We find the 
lack of silicate absorption and small CO column density toward DoAr 24Eb suggest the mid-infrared 
continuum is not as extinguished as the near-infrared, possibly due to the mid-infrared originating 
from an extended region. This, along with the velocity profile of the CO absorption, suggests the 
source of high extinction is likely due to a disk or disk wind associated with DoAr 24Eb. 
Subject headings: circumstellar matter - protoplanetary disks - stars: pre-main-sequence - stars: 
individual (DoAr 24E) 



1. INTRODUCTION 

Infrared Companions (IRCs; Koresko et al. 1997) are 
stellar objects characterized as being in binary systems 
with young stars, being faint in the visible, and dominat- 
ing the system flux at infrared and longer wavelengths. 
Being companions to young T Tauri stars suggests they 
are young stars as well, but the source of high extinction 
has remained a mystery. Koresko et al. (1997) noted the 
seven classical IRCs because they are well-studied sys- 
tems, so it is important to identify the sources of extinc- 
tion for future studies. As an example, if the extinction 
is due to IRCs having disks seen edge-on, these systems 
would present a good opportunity to search for molecular 
absorption originating in the inner disk regions. GV Tau 
N, a classical IRC, was found to exhibit warm molecu- 
lar C2H2 and HCN absorption (Doppmann ct al. 2008; 



cf. Gibb et al. 2007), rare features previously found in 
only one other disk system (IRS 46; Lahuis et al. 2006), 
which provided an important benchmark for inner disk 
chemistry models (e.g. Agundez et al. 2008; Woods & 
Willacy 2009). Another disk scenario that may cause 
the extinction seen toward IRCs is if the disk around the 
primary extends into the line of sight to the IRC (e.g. 
VV CrA; Smith et al. 2009, Kruger et al. 2011). Molec- 
ular absorption seen in this geometry would originate in 
the outer regions of the primary disk, likely resulting in 
largely different temperatures and column densities than 
expected toward an edge-on disk. Thus it would be im- 
portant to discern between these disk geometries before 
further studies of edge-on disks are confused by these 
disk systems. 

DoAr 24Eb is a classical infrared companion to the 
young star DoAr 24Ea, seen with a separation of 2". 03 
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with the IRC at a position angle of 150° (E of N) relative 
to the primary (Ghez et al. 1993). Koresko et al. (1997) 
modeled the primary as being a K0 type star with tem- 
perature 5240 K, while the IRC was modeled as having 
a central stellar temperature of 4850 K. We observed 
this system as part of our Spitzer Infrared Spectrograph 
(IRS) observing program of edge-on circumstellar disks 
and IRCs. We present here findings for the DoAr 24E 
system and investigate the disk geometry of this sys- 
tem and the reason for the high extinction toward DoAr 
24Eb. In section 2, we introduce the observations taken. 
In section 3, we describe the data reduction and discuss 
the detections. We then discuss modeling and conclu- 
sions in section 4, and give a summary of our findings in 
section 5. 

2. OBSERVATIONS 
2.1. Spitzer IRS Spectrum 

DoAr 24E was part of our Spitzer Cycle 5 General Ob- 
servers campaign, Program ID 50152, to use the Infrared 
Spectrograph (IRS; Houck et al. 2004) in short-low and 
short-high (SL and SH) stare mode. The SL module 
has a slitwidth of 3.6", plate scale of 1.8" pixel -1 (270 
AU pixel -1 assuming a distance of 150 pc to the system), 
and covers the spectral range 5.2-14.5 y.m with resolution 
R=60-127, while the SH module has a slitwidth of 4.7", 
plate scale of 2.3" pixel -1 (350 AU pixel -1 ), and covers 
9.9-19.6 ^m with R=600. The observations were designed 
to search for weak absorption features and so followed the 
observing technique described in Carr & Najita (2008). 
We shortened the exposure time of each data frame to 6 
seconds and increased the number of integrations to ver- 
ify repeatability of the observations. Each setting (SL1, 
SL2, and SH) had twelve on-source data frames at each 
nod position, accompanied with half as many off-source 
integrations of the blank sky for hot pixel subtraction 
and to remove background emission. The observation 
(AOR: 25679616) was taken on 2008 September 11 at 
16 ft 26 m 23.36 s - 24°21'1".6. The coordinates were cho- 
sen to center the slit on DoAr 24Eb, calculated by using 
the position of DoAr 24Ea from SIMBAD and the offsets 
found by Chelli et al. (1988). The observation used the 
nominal pointing accuracy of 0".l. 

The slit position angles during the observations were 
139.4° for the SH module and 2.6° for the SL. Given the 
primary-IRC position angle of 150° and separation 2". 03 
(Ghez et al. 1993), and the widths of the SH and SL slits 
(4". 7 and 3". 6, respectively), we would expect a major- 
ity of the flux from both binary components to be within 
the slit for both observations. A 4% continuum flux dif- 
ference between the SH and SL modules was corrected 
by scaling the SL spectrum by a constant, without any 
adjustments to the continuum shape, to match the SH 
spectrum at wavelengths where the orders overlap. 

2.2. Phoenix Spectrum 

On 2009 April 6, we took long-slit spectra with 
Phoenix (Hinkle et al. 2002) on Gemini South to exam- 
ine the CO fundamental vibrational transitions. We used 
the standard four pixel slit (0.35") with R=50,000, along 
with the M2150 filter with the echelle angle oriented to 
detect the 12 CO v=l-0 R(3)-R(5) transitions near 4.62 
[im, and the M2030 filter for the P(20)-P(21) transitions 



TABLE 1 

Fluxes as Measured with T-Recs Imaging 



Filter 


DoAr 24Ea (Jy) 


DoAr 24Eb (Jy) 


Si-1, 7.73 urn 


0.79±0.2 


1.96±0.3 


Si-2, 8.74 (im 


0.71±0.05 


1.94±0.06 


Si-4, 10.38 iim 


0.66±0.05 


1.92±0.07 


Si-6, 12.33 iim 


0.69±0.08 


1.88±0.1 


Qa, 18.3 [im 


1.16±0.5 


2.25±0.7 



near 4.86 uxn. Each ABBA nod position had a 120 sec- 
ond integration with four 30 second exposures co-added 
internally. Standard star observations were taken before 
the science frames for telluric division. Only DoAr 24Eb 
was within the slit during the observation due to a sign 
error in the position angle given for the queue observa- 
tion. 

2.3. T-ReCS Imaging 

To search for extended emission and to find flux ratios 
between binary components, we used T-ReCS (Telesco et 
al. 1998) on Gemini South to image DoAr 24E in 6 filters 
ranging from 7.7-18.3 [im (see Table 1). We took calibra- 
tion images of standard stars before and after DoAr 24E. 
Observations were taken by cycling through all six filters 
once on target. We used the standard 3 position chop- 
nod observing sequence and used the off-beams only for 
background subtraction. These observations were taken 
2009 June 8. 

3. DATA REDUCTION AND DETECTIONS 

3.1. Spitzer IRS Spectrum 

For both the high resolution SH and low resolution SL 
spectra, we created custom IDL routines to identify and 
remove bad pixels on the array found using the images of 
the nearby blank sky. The SL spectra were then reduced 
using optimal point source extraction in the Spitzer IRS 
Custom Extractor (SPICE; version 2.2). The high reso- 
lution SH spectra were reduced following the reduction 
technique described in Carr & Najita (2008) using cus- 
tom IDL routines along with standard spectral reduction 
routines in IRAF. We first used the apall routine in IRAF 
to extract the raw spectra. We used the standard star 
HR 6688 for flux calibration and to remove fringing. Be- 
cause fringing is dependent on the target position in the 
slit, we used the observation of HR 6688 that resulted in 
the least amount of fringing at shorter wavelengths after 
division. We then used IRSFRINGE in IDL to remove 
any residual fringing. The final spectrum is shown in 
Figure 1. 

The Spitzer IRS spectrum shows warm molecular H 2 0, 
C2H2, HCN, and CO2 emission as was found by Pontopp- 
idan et al. (2010). We examined the archival Spitzer IRS 
data reported by Pontoppidan et al. (2010) and do not 
detect any significant changes in emission. Further emis- 
sion studies of DoAr 24E will be carried out with ground- 
based data by Moerchen et al. (in prep.). The system also 
shows a double-peaked spectrum in the 8-12 [im region. 
The shape of this spectral region is similar to SVS20 (cf. 
Figure 7; Alexander et al. 2003), which was found to be 
due to hot, optically thin silicate emission being extin- 
guished by cool, optically thin silicate absorption. Such a 
combination of silicate emission and absorption is a com- 
mon interpretation for the relatively flat spectra found 
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Fig. 1. — Spitzer IRS spectrum of DoAr 24E, combining both binary components. Overplotted are the combined photometric fluxes 
measured from our T-ReCS imaging, summing 95% of the flux from the DoAr 24E components. The inset indicates the C2H2, HCN, and 
CO2 emission features. 



in YSOs (Mitchell & Robinson 1981; Furlan et al. 2008). 
Alternatively, the peaks at 9.2 and 11.3 |im may be due 
to crystalline silicates which are also commonly seen in 
YSOs (Olofsson et al. 2009). The lack of a strong silicate 
absorption feature is surprising given the previous high 
extinction estimates (i.e. 26 mag; Koresko et al. 1997). 
Using the relation Ay/rsio = 16.6 ± 2.1 from Rieke & 
Lebofsky (1985) for grains in the ISM, we would expect 
to see a silicate absorption feature with tsw = 1.6 ± 0.2. 

3.2. Phoenix Spectrum 

The near-IR Phoenix spectra of DoAr 24Eb were re- 
duced in IRAF. We first used the identify task on the sky 
emission lines to get a wavelength solution and to rectify 
the slit in each frame, with the rest frame wavelengths 
of the emission lines taken from the HITRAN database 
(Rothman et al. 1998). 

In the Phoenix spectrum, we found 12 CO transitions 
both in emission and absorption, shown in Figure 2. The 
emission component at least partially arises from warm 
gas as we see the v=2-l R(ll)-R(13) and P(15) transi- 
tions, as well as the high-J v=l-0 P(21) transition. There 
are two discrete absorption features, found exclusively 
in the low-J R(3)-R(5) transitions, at vlsr = —5.8 km 
s" 1 (hereafter A-6) and 3.2 km s" 1 (A3). We fit the 
normalized, average profile of the R(3)-R(5) transitions 
using Gaussian distributions for both the emission and 
absorption features. The fit is shown in Figure 3, and 
the Gaussian parameters used are given in Table 2. In 
calculating the equivalent widths, given in Table 3, we 
assume the relevant local continuum includes the emis- 
sion component as appropriate if the emission is coming 
from the innermost regions of the disk and the absorbing 
gas is further out. We also assume the covering factor, 
the fraction of the unresolved emitting area that is being 



TABLE 2 

Measured Parameters for the 12 CO Line Profiles 







visa (km s -1 ) 


FWHM (km s- 1 ) 


Emission 




6.1 ±0.6 


76 ±2 


Absorption 


(A-6) 


-5.8 ±0.2 


6.8 ± 0.3 


Absorption 


(A3) 


3.2 ±0.3 


6.1 ± 0.8 



TABLE 3 

Measured Equivalent Widths of the 12 CO Fundamental 
Absorption Lines 

Component Line ID Equivalent Width (cm -1 ) 



A3 



A-6 



R(3) 
R(4) 
R(5) 
R(3) 
R(4) 
R(5) 



0.0173 ±0.001 
0.0178 ±0.001 
0.0160 ±0.001 
0.0145 ± 0.0009 
0.0155 ± 0.0009 
0.0143 ± 0.0009 



absorbed, to be 1. 

3.3. T-ReCS Imaging 

We used custom IDL routines to reduce the T-ReCS 
images. High winds during the observation greatly re- 
duced the image quality, so the images were centered 
on the peak of the target flux before combining. We 
searched for extended emission by varying the extraction 
apertures, but found none. Deconvolution using lucy in 
IRAF gave the same results. 

We detected both binary components in all images. 
We used a 1" aperture radius to estimate the photomet- 
ric fluxes for the individual binary components, and the 
results are listed in Table 1. We find a binary separa- 
tion of 2".03 ± 0".01 at a position angle of 149°.0 ± 0°.l. 
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Fig. 2. — CO fundamental spectra of DoAr 24Eb as observed by Phoenix in the topocentric frame in orders 16 (above) and 15 (below) 
with detected transitions marked. The telluric transmission spectrum from a standard star is also plotted below each science spectrum. 
The dotted line in order 15 is the average velocity profile fit to the R(3)-R(5) transitions, shown in Figure 3, centered at the P(21) transition 
for comparison. 
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Fig. 3. — The normalized, average CO line profile from the 
Phoenix using the R(3)-R(5) transitions. The solid line is a fit 
using Gaussian profiles with the parameters listed in Table 2. The 
missing data points are due to telluric absorption. 



Ghez et al. (1993) also measured a separation of 2". 03 in 
July 1990, so we conclude that these objects have similar 
proper motions and are associated with each other. As 
the change in separation since this observation is < 0".01, 
this shows the velocity difference between the objects in 
the direction perpendicular to our line of sight is <3.7 
km s _1 , assuming the system is 150 pc away. We also 
find the combined fluxes are similar to the Spitzer IRS 
spectrum, shown in Figure 1, assuming a 5% flux loss 
from both components during the IRS observation. We 



assume a flux loss from both components as the IRS slit 
was within 10° of the actual position angle of the DoAr 
24E components during the observation. 

4. DISCUSSION 

4.1. CO Gas 

Najita et al. (2003) showed that CO emission is com- 
monly detected toward T Tauri stars, and most likely 
originates from the inner 5 AU of the disk. They found 
the emission lines typically had FWHM of « 70 km 
s _1 , similar to that found in the R(3)-R(5) transitions of 
DoAr 24Eb (FWHM= 76 ±2 km s" 1 ). We thus continue 
with the assumption that this emission largely originates 
from the inner 5 AU of the disk of DoAr 24Eb. We do 
note, however, the P(21) transition has a different ve- 
locity profile than the R(3)-R(5) transitions as it peaks 
at -15 km s _1 , has a FWHM of 45 km s _1 , and is not 
Gaussian. While the v=2-l P(15) transition complicates 
the profile, the red side of the v=l-0 P(21) transition 
appears to be broader than the blue side. 

We modeled the CO gas absorption components to ex- 
tract temperature and column density estimates. The 
velocity widths of these components are near the best- 
case instrument resolution of 5.88 km s _1 (Hinkle 1999), 
so we cannot get precise measurements, but we can find a 
range of values that would result in the observed equiva- 
lent widths. We modeled the transitions using a uniform 
temperature slab model with a covering factor of 1 and 
assuming the gas is in local thermodynamic equilibrium. 
We considered a wide range of temperatures and column 
densities and used a x 2 goodness-of-fit to find the best-fit 
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Fig. 4. — The goodness-of-fit contours for the temperature and 
column density for the two absorption components, indicated by 
the velocity shifts of the corresponding components. The contours 
show the 68% and 90% confidence levels assuming the best-case 
instrument resolution of 5.88 km s _1 (see text). 

parameters. The results are shown in Figure 4. 

The CO absorption feature at vlsr — 3.2 km s -1 (A3) 
has a small FWHM of 6.1 km s _1 . Using the best-case 
instrument resolution of 5.88 km s -1 , we modeled the 
absorbing gas with an intrinsic velocity width of 1.6 km 
s -1 . The temperature and column density are not well 
constrained; the temperature we found is 75±50 K with 
a column density of 5.5 ± 4.0 x 10 17 cm~ 2 , but these 
parameters are sensitive to the turbulent velocity. If we 
assume an instrument resolution of 6.0 km s^ 1 , just less 
than the velocity width of A3 and corresponding to an 
intrinsic turbulent velocity of 1.1 km s _1 , the estimated 
l-cr column density increases to 2.5 i 1.8 x 10 19 cm -1 
and the temperature decreases to 33 ± 10 K. 

The absorption feature at vlsr = —5.8 km s _1 (A-6) 
has a FWHM of 6.8 km s _1 . We modeled these transi- 
tions with a turbulent velocity of 3.4 km s" 1 and found 
a best- fit temperature of 115 ± 30 K and column den- 
sity 5.5 ± 0.3 x 10 16 cm~ 2 . These estimates show little 
dependence on the turbulent velocity used. Again us- 
ing an instrument resolution of 6.0 km s _1 and resulting 
turbulent velocity of 3.2 km s" 1 , we find a temperature 
estimate of 105±20 K and column density 5.9±0.3 x 10 16 
cm~ 2 . 

To further investigate the CO absorption, we ob- 
tained spectroscopic data from the CRyogenic high- 
resolution InfraRed Echelle Spectrograph (CRIRES; 
Kaufl et al. 2004) that was reported in Pontoppidan et 
al. (2011). These spectra include the 12 CO R(0)-R(1) 
and P(l)-P(9) transitions and were observed on 2007 
September 3. The spectrum of DoAr 24Ea shows only 
the A3 absorption lines. The emission line profile for 
DoAr 24Eb does not appear to change between the low- 
and high-J transitions, up to P(9). We also found the 
A3 transitions have an average FWHM width of 3.6 km 
s _1 , while the A-6 lines have an average width of 4.9 km 
s" 1 . Given the measured instrument resolution of 3.18 
km s" 1 (CRIRES User Manual, Issue 85.2), this would 
indicate FWHM turbulent velocity widths of 1.7 and 3.7 
km s _1 , respectively. This is comparable to the turbu- 
lent velocity estimates with the Phoenix data (1.6 and 
3.4 km s" 1 , respectively). The A3 line is only seen up 
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Fio. 5. — The normalized, average CO line profile from the 
CRIRES spectrum using the P(3)-P(5) transitions. The solid line 
is a fit to the data with the emission estimated as a Gaussian with 
FWHM=54 km s _1 and vlsr = —9 km s _1 . The two absorp- 
tion lines were found by using the turbulent velocity estimates and 
average absorption strengths found in the Phoenix spectrum, con- 
volved with the CRIRES resolving power of 3.18 km s (see text), 
showing a good fit to the absorption lines in the CRIRES spectrum. 

to the P(7) transition, showing it originates in cooler gas 
(<50 K). It is difficult to estimate the optical depth of the 
A-6 transitions due to being so near the strong telluric 
CO absorption and the emission line profile changing for 
higher transitions. However, the A-6 absorption lines are 
seen through the P(9) transition, and we find that mod- 
eling the gas with temperature 115 K and column density 
5.5 x 10 16 cm -2 , as found with the Phoenix spectrum, 
fits the CRIRES spectrum well. 

The emission line in the CRIRES spectrum is slightly 
narrower than found with Phoenix, and peaks at nega- 
tive velocities. If we again estimate the emission as a 
Gaussian, we find it has a FWHM of 54 km s _1 and 
vlsr = —9 km s" 1 (see Figure 5). This change in veloc- 
ity from the emission seen in the Phoenix spectrum may 
be due to complicated details within the inner regions of 
the disk, possibly due to DoAr 2Eb itself being a binary 
(Koresko 2002). 

4.2. Gas Location 

There are several possible locations for the two absorb- 
ing CO gas columns, among which are a disk, disk wind, 
or envelope associated with DoAr 24Eb, the outer regions 
of a disk around DoAr 24Ea, or the interstellar medium 
(ISM). 

The main source of interstellar extinction is likely the 
p Ophiuchus (p Oph) cloud. This cloud has a system 
velocity of vlsr. ~ 3 — 4 km s _1 , similar to the A3 CO. 
Chelli et al. (1988) estimated an interstellar extinction of 
Ay =6.2 mag for DoAr 24Ea, so if we use the interstel- 
lar ratio of N(CO) /Ay = 1.4 x 10 17 cm" 2 mag" 1 (Rettig 
et al. 2006), we would expect a CO column density of 
8.7 x 10 17 cm -1 . This column density is within the 1-a 
error of the estimate for the A3 CO, and the low temper- 
ature of the absorbing gas (<50 K) is typical for the ISM. 
The CRIRES spectrum also includes DoAr 24Ea within 
the slit (although the spectrum of DoAr 24Ea was not 
included in Pontoppidan et al. 2011) and we find that 
DoAr 24Ea also shows the A3 line with a similar velocity 
profile and equivalent width as DoAr 24Eb. We conclude 
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Fig. 6. — The best-fit models to the SEDs of DoAr 24Ea and DoAr 24Eb. The dots are spectral points taken from our T-ReCS imaging 
and Chelli et al. (1988). The solid lines are the SED fits, while the dashed lines are the stellar photospheres as it would look without a 
circumstellar disk but still including the interstellar extinction. Left: The spectrum of DoAr 24Ea. Middle: The spectrum of DoAr 24Eb 
modeled as having an interstellar extinction <10 mag. Right: The spectrum of DoAr 24Eb modeled as having an interstellar extinction 
<30 mag. 



the A3 CO gas is associated with the ISM. 

The A-6 line is not seen toward DoAr 24Ea in the 
CRIRES spectrum, so this CO is associated with the 
material that is extinguishing DoAr 24Eb but not DoAr 
24Ea. Koresko et al. (1997) used H and K band photom- 
etry to model the infrared companion DoAr 24Eb and 
estimated a high extinction of Ay = 26 mag. Removing 
the interstellar extinction, this would indicate Ay « 20 
mag for the material that is obscuring DoAr 24Eb but 
not DoAr 24Ea. Again comparing to the N(CO)/Ay ra- 
tio of the ISM, this would indicate a high CO column 
density of 2.8 x 10 18 cm -2 , which is nearly 2 orders of 
magnitude greater than what we estimated for the A-6 
CO. While this implies a substantially smaller gas/dust 
ratio than is seen in the ISM (« 1.5%), it may be because 
the M-band continuum is coming from a larger region of 
the disk than where the near-infrared extinction is oc- 
curring, resulting in a covering factor of <1. This may 
also explain why we do not see a large silicate absorp- 
tion feature in the Spitzer IRS spectrum as expected for 
a system showing such high extinction. 

While the location of the A-6 CO is unclear, it would 
not be associated with the outer disk regions of DoAr 
24Ea. Assuming a distance of 140-160 pc to the sys- 
tem, the binary components have a projected separation 
of 285-325 AU, and the gas temperature of ~ 100 K is 
high for these disk radii. Furthermore, Doppmann et 
al. (2003) observed the photospheric absorption of DoAr 
24Ea (GY 20A) to find a redshift of v LSR = 6.4 km s" 1 . 
If the A-6 CO was associated with DoAr 24Ea, the gas 
would have a 12.2 km s _1 velocity shift from the central 
stellar component, much too high for gas at these disk 
radii. With such a large column only seen toward DoAr 
24Eb, and having a discrete velocity shift of -9 km s _1 
from the surrounding cloud, it is difficult to envision a 
physical situation in which this is not associated with the 
disk of DoAr 24Eb. It is unknown if the A-6 line is at the 
system velocity of DoAr 24Eb, or possibly shifted from 
the central star velocity, such as in a disk wind. 

4.3. SED Modeling 

To further investigate the system, we used the spec- 
tral energy distribution (SED) model fitter described in 
Robitaille et al. (2006). This tool finds a best-fit to an 
input SED using a database of 200,000 pre-computed ra- 
diation transfer models of circumstellar disks, with 14 



physical parameters varied at 10 different disk viewing 
angles ranging from 18.2 to 87.1°. The models differenti- 
ate the circumstellar extinction, due to the disk material 
between the stellar surface and the outer edge of the disk, 
and the interstellar extinction between the disk edge and 
the observer. The user inputs spectral points, and can 
constrain the distance to the system and the interstellar 
extinction. The best-fit models are then returned in or- 
der of increasing \ 2 value. The estimates for the distance 
to the p Oph star forming region are in the 140-160 pc 
range, we found that using this limited distance range re- 
sulted in only a few returned models. We thus extended 
the distance range by 20 pc in both directions (120-180 
pc) in order to get a range of disk parameters. We then 
used only those models with X^ggj- — X 2 < 5N, where N 
is the number of input data points, and stellar tempera- 
tures within 500 K of estimates in Koresko et al. (1997). 
We show the best-fit models described below in Figure 
6. Table 4 reports parameters for the best-fit models, as 
well as the minimum and maximum values of each indi- 
vidual parameter for the models returned given the \ 2 
cut-off. 

We first fit the SED of DoAr 24Ea for comparison with 
predictions. We constrained the interstellar extinction to 
be in the 0-10 mag range, and find the models are close to 
previous estimates. While Koresko et al. (1997) modeled 
this system as having a central stellar temperature of 
5240 K, the model stellar temperatures are consistently 
lower (« 4800 K). This may account for the modeled 
interstellar extinction (5.3 mag) being slightly lower than 
previous estimates (e.g. 6.2 mag; Chelli et al. 1988). 
The models predict the system does not have an infalling 
envelope. 

To investigate the scenario where the extinction in 
DoAr 24Eb is due to a circumstellar disk, likely seen near 
edge-on, we again constrained the interstellar extinction 
to be in the 0-10 mag range. The returned models have 
interstellar extinctions that are higher than predictions, 
tending to be at the upper limit of 10 mag. The best-fit 
model also surprisingly predicts a small circumstellar ex- 
tinction Ay of only 0.026 mag. The extinction estimates 
go as high as 150 mag, but the median extinction esti- 
mate in the returned models is 12 mag, closer to expec- 
tations given interstellar extinction constraint and pre- 
vious total extinction estimates (i.e. 26 mag; Koresko et 
al. 1997). Unlike DoAr 24Ea, most of these models also 



Gas and Dust Absorption in the DoAr 24E System 



7 



TABLE 4 

Estimated Disk Parameters from SED Fitting 







DoAr 24Eb 


DoAr 24Eb 


Parameter 


DoAr 24Ea 


(Ay =0-10 mag) a 


(Ay = 0-30 mag) a 


Interstellar Ay (mag) 


5.3 (5.2 5.3) 


10 (1.6 10) 


17 (11 24) 


Circumstellar Ay (mag) 


2.2e-4 (2.2e-4 6.8e-4) 


0.026 (0.026 ISO) 6 


2.5 (0.0007 4.6) 


Stellar Temperature (K) 


4810 (4710 4810) 


4550 (4400 5100) 


4600 (4400 5100) 


Stellar Mass (Mq) 


2.0 (1.8 2.0) 


1.5 (1.3 3.1) 


1.9 (1.5 3.5) 


Env. Infall Rate (1CT 6 M yr" 1 ) 


(0 0) 


5.5e-3 (0 13) 


2.1 (0 3.8) 


Disk Accretion Rate (10~ 8 M yr" 1 ) 


8.3 (8.3 11) 


86 (0.26 86) 


6.1 (6.3e-3 86) 


Disk Inner Radius (AU) 


0.15 (0.15 0.15) 


1.4 (0.21 1.6) 


0.63 (0.26 2.3) 


Cavity Opening Angle ( ° ) 


No env. 


32 (19 40) c 


35 (32 54) c 


Disk Inclination ( ° ) 


76 (63 81) 


76 (63 81) 


70 (32 81) 



Notes. Shown are the best-fit and the (min, max) of the returned parameters for the best-fit models. 

"Interstellar extinction constraint. 

6 Median value of 12 mag, see text. 

c The range for those disks with an envelope. 



predict an infalling envelope. Given that DoAr 24Eb is 
physically associated with DoAr 24Ea, as shown by the 
common proper motions, we would expect the systems 
to be the same age and likely the same stage of evolu- 
tion. Thus, having an envelope around DoAr 24Eb and 
not DoAr 24Eb would require some mechanism to either 
prolong the envelope around DoAr 24Eb or shorten the 
envelope around DoAr 24Ea. On a cautionary note, the 
SED model fitter only tests disk inclinations at discrete 
intervals, which is a drawback for edge-on disks where 
small inclination changes result in large SED changes 
(Robitaille et al. 2007; Kruger et al. 2011). Presumably, 
finer gradations in disk inclination would provide models 
with more appropriate circumstellar extinction estimates 
and possibly eliminate the envelopes. 

For completeness, we modeled DoAr 24Eb as having 
another compact, interstellar component in the line of 
sight, as would be the case if the outer disk of DoAr 
24Ea was in the line of sight, by allowing the interstel- 
lar extinction to be in the range 0-30 mag. Because 
there were so many best-fit models returned, we lowered 
the x 2 cut-off to x? , — x 2 < 3N so that we were not 
over-estimating our parameters (Robitaille et al. 2007). 
The stellar temperatures are comparable to Koresko et 
al. (1997), ranging above and below, and most models do 
not include a surrounding envelope. The modeled inter- 
stellar extinctions are slightly lower but similar to pre- 
vious extinction estimates for DoAr 24Eb (i.e. 26 mag; 
Koresko et al. 1997). The models associated a majority 
of the extinction with the ISM, with the largest circum- 
stellar extinction being 0.44 mag. As discussed in section 
4.2, it is unlikely that such a large column can be seen 
toward DoAr 24Eb, and not DoAr 24Ea, without being 
a part of the disk of DoAr 24Eb. 

5. SUMMARY 

We investigated the source of the high extinction to- 
ward the infrared companion DoAr 24Eb, and considered 
several disk and envelope geometries that may explain 
the high extinction toward DoAr 24Eb: a disk or en- 
velope around DoAr 24Eb, the outer regions of a disk 
around DoAr 24Ea, or the interstellar medium (ISM). 
To test the scenarios, we analyzed near-infrared spec- 
tra from Phoenix of the CO fundamental transitions and 
modeled the mid-infrared photometric fluxes we found 



with T-ReCS imaging. 

DoAr 24Eb shows two 12 CO absorption features at 
vlsr = 3.2 (A3) and -5.8 km s" 1 (A-6). Using a slab 
model and optically thick gas, we found the A-6 CO has 
a temperature of ~ 100 K and column density ~ 6 x 10 16 
cm~ 2 while the A3 CO has a higher column density 
(10 17 - 10 18 cm- 2 ). A CRIRES spectrum also indicates 
the A3 line has a low temperature (<50 K). The velocity 
shift, column density, and temperature of the A3 line, as 
well as its presence in the DoAr 24Ea spectrum, indicate 
that it is associated with the interstellar medium. 

Since A-6 is seen toward DoAr 24Eb and not DoAr 
24Ea, it is most likely associated with DoAr 24Eb. This 
CO is blueshifted ~ 12 km s _1 from the radial velocity of 
DoAr 24Ea (Doppmann et al. 2003) and has a tempera- 
ture of <~ 100 K, so it would not be in outer disk regions 
of DoAr 24Ea in the line of sight to DoAr 24Eb. The 
discrete shift from the cloud velocity and having such a 
large column in the line of sight toward only DoAr 24Eb 
make it most likely in a disk or disk wind associated with 
DoAr 24Eb. 

The Spitzer IRS spectrum does not show strong sili- 
cate absorption as we would expect for an object showing 
high extinction. Also, our measurement of the column 
density for the A-6 CO absorption, which is associated 
with DoAr 24Eb, is two orders of magnitude smaller than 
expected for the high extinction. These factors suggest 
the mid-infrared continuum is not as extinguished as the 
near-infrared. This may be due to the source of the mid- 
infrared flux being more extended than the near-infrared, 
possibly due to an extended disk. This further argues 
against the source of extinction toward the IRC being 
due to the disk around the primary or the ISM, which 
would be expected to show a more uniform extinction. 

To further investigate the disk scenarios, we modeled 
the spectral energy distributions of both components. 
We tested if the extinguishing material is likely associ- 
ated with the disk around DoAr 24Eb by constraining the 
interstellar medium to be less that 10 mag, as is found for 
DoAr 24Ea. We found the models predict a wide range 
of circumstellar extinctions and an infalling envelope for 
DoAr 24Eb, which is not predicted for DoAr 24Ea. How- 
ever, this may be due to the models only being tested at 
discrete intervals (Robitaille et al. 2007), which results 
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in large SED changes for high inclinations. 
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